First stereoscopic coronal loop reconstructions from 
STEREO/SECCHI images 
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ABSTRACT 

We present the first reconstruction of the three-dimensional shape of magnetic loops in an 
active region from two different vantage points based on simultaneously recorded images. The 
images were taken by the two EUVI telescopes of the SECCHI instrument onboard the recently 
launched STEREO spacecraft when the heliocentric separation of the two space probes was 12 
degrees. We demostrate that these data allow to obtain a reliable three-dimensional reconstruc- 
tion of sufficiently bright loops. The result is compared with field lines derived from a coronal 
magnetic field model extrapolated from a photospheric magnetogram recorded nearly simultane- 
ously by SOHO/MDI. We attribute discrepancies between reconstructed loops and extrapolated 
field lines to the inadequacy of the linear force-free field model used for the extrapolation. 



Subject headings: solar corona, magnetic field, stereoscopy 



1. Introduction 

With the launch of NASA's STEREO mission 
in October 2006, a new dimension of solar coronal 
observations has been opened. For the first time, 
objects above the solar surface can be perceived 
in three dimensions by analysing the stereo image 
pairs observed with the SECCHI instruments on- 
board the STEREO spacecraft and without mak- 
ing a-priori assumptions about their shape. The 
two STEREO spacecraft orbit the Sun at approx- 
imately 1 AU near the ecliptic plane with a slowly 
increasing angle of about 45 degrees/year between 
STEREO A and STEREO B. Each spacecraft is 
equipped with, among other instruments, an EUV 
telescope (SECCHI/EUVI). For the objectives of 



the mission and more details about the EUVI tele- 
scopes see 



(2007) 



Wuelser et al. (2004) and Howard et al 
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The major building blocks of the solar corona 
are loops of magnetic flux which are outlined by 
emissions at, e.g., EUV wavelengths. In princi- 
ple, the magnetic field in the lower corona can 
be derived from surface magnetograms by way of 
extrapolations (e.g. |Wiegelmanii]|2007| . However, 
missing boundary values and measurement errors 
may introduce considerable uncertainties in the 
extrapolation results so that there is an obvious 
need for an alternative three-dimensional deter- 
mination of the coronal magnetic field geometry. 
Among other goals of the mission, this require- 
ment has been one of the drivers for STEREO. 

Attempts for a three-dimensional reconstruc- 
tion of the coronal magnetic field from EUV ob- 
servations have started long before STEREO data 
was available and date back more than a decade 
( Berton fc Sakurai|1985 Kouchmy & Molodensky 



1992). Here, we for the first time use two simul- 



taneously observed EUVI images observed by the 
two STEREO probes and rigourously reconstruct 
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loop shapes without any further assumption about 
their temporal or spatial behaviour from which 
earlier reconstructions employing consecutive im- 



ages from a single spacecraft suffered (Aschwan- 



den et al. 2006 1 . We compare the reconstruction 
results with field lines derived from linear force- 
free magnetic field models with variable a, the ra- 
tio of field-aligned current density to field strength 
QSeehafer|[l9#l ). 



2. The data 

For our reconstruction we used EUV images 
at A = 17.1 nm taken by the almost identical 
SECCHI/EUVI telescopes onboard of the two 
STEREO spacecraft at 2007-06-08 03:21 UT when 
the well isolated active region NOAA 0960 was 
close to solar disk centre. The line A = 17.1 nm is 
emitted by the Fe IX ion which in thermal equilib- 
rium forms at about 1.1 million K. At the time of 
these observations, the two STEREO spacecraft 
had a heliocentric separation of 11.807 degrees. 
The precise spacecraft positions at the time of the 
observation are listed in table Q] 

For a comparison of our reconstruction with 
magnetic field lines we made use of a SOHO/MDI 
magnetogram (Scherrer et al. 1995[ ) taken only 
9 seconds prior to the EUVI images. The ac- 
tive region is well isolated from neighbouring field 
sources so that an extrapolation of the surface field 
is possible. MDI, however, provides only the line- 
of- sight field component, which for this bipolar re- 
gion close to the disk centre is almost identical to 
the radial field component on the solar surface. 
For this reason we can employ here only a linear 
force-free field model for the extrapolation of the 
magnetogram ( |Seehafer|[T978| . 



STEREO probe 



B 



Helioc. dist. (AU) 
Sun's app. rad. (arcsec) 
Longitude (degrees) 
Latitude (degrees) 



1.068788 0.958071 

897.866 1001.625 

-4.277 7.524 

-0.293 0.095 



Table 1: STEREO spacecraft coordinates at the 
time of the observations. Spacecraft longitude and 
latitude are given in the Heliocentric Earth Eclip- 
tic (HEE) coordinate system. 



3. The reconstruction 

The first step in the stereoscopic reconstruc- 
tion scheme is the isolation and identification of 
individual loops in each of the EUV images. In 
Figure [l] we show the portion of the EUV images 
containing the active region. The EUV structures 
were contrast enhanced by an unsharp mask filter. 
Next, individual loop structures were detected by 
a loop segmentation program. This program de- 
tects individual bright loops in an image by treat- 



ing them as elongated intensity ridges (Inhester 



et al. 20071. For identification, the loop curves 
were enumerated. These assignments, e.g. a num- 
ber ia for a loop curve in image A, were made 
independently in each image. 

To establish correspondences of projections 
*A <^ t-b °f the same loop across the images is 
the hardest part in the stereoscopy procedure. 
For isolated loops they can sometimes be guessed 
by visual comparison of the image pair. Also, 
some guidance is provided by matching constraints 
which corresponding pairs of loop projections have 



to obey (Inhester 2006). Often, however, the vi- 
sual comparison of loop structures does not yield 
unique correspondences. To disentangle the typ- 
ically crowded active region loop ensembles we 
have developed a systematic scheme which deter- 
mines correspondences with the help of magnetic 



field model calculations 


Wiegelmann et al. 2005 


Wiegelmann & Inhester 


2006 




Feng et al. 


2007). 



The idea is to find three dimensional field lines 
from a more or less accurate model of the ac- 
tive region magnetic field as a first approximation 
to the final loops whose projections are close to 
the loop projections identified in the images from 
spacecraft A and B. If a field line can be found 
with projections sufficiently close to a loop in both 
images, this is strong evidence that these loop 
curves represent projections of the same three- 
dimensional loop. 

We quantify the proximity of a projected field 
line / to a loop curve %a in image A, say, by 
the mean distance Ca(ia,1) between the two- 
dimensional curves in this image. The probabil- 
ity of a correspondence between a pair {ia^b) of 
loop curves in image A and B can then be mea- 
sured by C = ^min; (Ca{ia, I) + Cs(is) 0)- H ere i 
the set of possible field lines I comprised all pos- 
sible foot point locations and a wide range of a 
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Fig. 1. — Contrast enhanced zoom of the EUVI images B (left) and A (right) of the active region NOAA 
0960. Heliographic north is upward. The axes are scaled according to the image pixel size. Individual loop 
structures are emphasized by white curves and enumerated. Equal numbers do not imply a correspondence 
across the images. 



1 /c (1/Mm) 




Fig. 2. — Proximity of loops identified in images A 
and B in figure [l] The proximity is expressed by 
the inverse of a distance measure C (see text). The 
loops from image A and B are arranged along axes 
'EUVI-A.' and 'EUVIJ3' according to their respec- 
tive identification number %a and is- For each pair 
(*-Ai*b) the inverse of C is displayed by a column 
at the location of the loop pair in this matrix rep- 
resentation. Columns exceeding 0.5 Mm -1 have a 
black top. Here, only loops from the northern half 
of the active region were considered. 



values from -0.01 to +0.01 Mm" 1 . The field lines 
I here only serve as a means to establish the cor- 
respondence, they are not intended to represent a 
consistent field model of the active region. The lin- 
ear force-free field model used is only consistent if 
a is a global constant. Strictly speaking, the field 
lines i m i n for which C attains the minimum are 
each from a different field line model as a turned 
out to differ for each loop pair. 

In Figure [2] the inverse of C is shown for the 
loops in the northern half of the active region, %a 
= to 18 and is = to 17. Some few loop combi- 
nations show a clearly enhanced 1/C and are thus 
much more probable than the majority of combi- 
nations {%Ai *b)- We accepted for a reconstruction 
only loop pairs with a value of C below 2 Mm. 
This corresponds to an average distance between 
the field line projection and the loop curves in each 
image of 2 pixels or less. When more than one 
combination was possible for one loop, the most 
probable one was taken such that each loop re- 
ceives no more than one partner and the sum of 
C of all selected correspondences was minimised 
(Wiegelmann & Inhester 2006). In all, 20 pairs 
from Figure [I] could thus be identified. 

The last step is the stereoscopic reconstruc- 
tion of the three-dimensional loop from each ac- 
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cepted pair This purely geometrical step 



often yields multiple solutions (Inhcster 2006) 



They were discarded by retaining only the three- 
dimensional reconstruction closest to the best fit 
field line Z m j n . 

4. Results 

In Figures [3] and [4] we present two views of a set 
of reconstructed loops. Figure [3] shows the recon- 
structed loops (yellow) and the associated closest 
fit field lines (red) obtained by extrapolation from 
a position within a degree from the STEREO A 
spacecraft. As expected, loops and field lines agree 
relatively well from this perspective because they 
were chosen to be close in this projection. Figure [4] 
therefore provides a completely different view of 
the active region. This view shows that most of 
the loops cannot easily be approximated by planar 
curve segments. This geometrical simplification 
was often used for loop reconstructions in the past 
because a more involved shape could only rarely 
be resolved from previous observations. This fig- 
ure also reveals deviations between the loops and 
field lines. E.g., the loops on presumably open 
field lines appear to be more strongly curved than 
the corresponding field lines from the extrapola- 
tion. 

We attribute this disagreement to a deficiency 
of the linear force-free field line extrapolation. For 
the closed field lines, the best fit \a\ values derived 
above fell in the range from 1.8 to 8.3 10~ 3 Mm -1 
(see table [2j . For the open field lines, these values 
turned out to be smaller in magnitude, with values 
\a\ <2.5 10 -3 Mm -1 . As a is a global constant 
for the linear force-free field model, the influence 
of the stronger currents on the closed active-region 
field lines is not accounted for on the open field 
lines. This may explain why the open field lines 
were calculated with less curvature than the cor- 
responding stereoscopically reconstructed loops. 

The loop reconstruction is also prone to errors, 
however. These may occur whenever a projected 
loop section in the images are directed tangentially 
to an epipolar line ( Inhester]|2006 1 . For the view- 
ing geometry of our observations, epipolar lines 
are nearly horizontal in the images and the crit- 
ical part for closed, E-W orientated loops there- 
fore lies more or less near their apex. Also the 
open loop structures 16-19 in image B and 17-20 




Fig. 3. — Vertical View of the three-dimensional 
reconstruction results from a viewpoint within a 
degree from the STEREO A spacecraft. Helio- 
graphic north is upward. The reconstructed loop 
sections are drawn in yellow, the closest fit field 
lines in red. The loop pairs {ia^b) drawn are: 
4-2, 12-12, 5-3, 7-5 (northward part of this AR) 
and 45-45, 44-43, 42-42, 24-23, 30-29 (southward 
part). 




Fig. 4. — Same as Figure [3] but seen from a view 
point NE of the active region. Heliographic north 
points to the lower left corner. The SECCHI in- 
struments observed from approximately above. 
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Fig. 5. — Example of a reconstructed loop with 
error estimates. The reconstruction is shown in 
yellow along with error bars. The associated best- 
fit linear force-free field line (red curve) is much 
lower in height. 

in image A (see Figure [TJ suffer from this problem 
as they arc orientated almost entirely horizontally 
in the images. We have therefore not attempted 
to reconstruct them even though a correspondence 
could well be identified. 

In Figure [5] we display the reconstruction of 
loop (5,3) (yellow curve) which shows by far the 
largest deviation to its best fit linear force-free 
field line (red curve). For most other loops, this 
discrepancy is much less although the agreement 
is rarely perfect. For some points along the loop 
(5,3), we also show error bars which represent the 
geometrical reconstruction error when the uncer- 
tainty for the loop projection in the images is as- 
sumed to be 1.5 pixels. In this case, the height of 



Loop pair 


\a\ 


height 


length 


M, is 


(lO^Afra" 1 ) 


(Mm) 


(Mm) 


5, 3 


1.8 


71.9 


229 


7, 5 


8.3 


20.6 


105 


45,45 


2.3 


58.2 


253 


44,43 


2.8 


27.3 


188 


42,42 


2.8 


57.2 


210 



Table 2: Best fit field line parameters for a rep- 
resentative list of closed loops of active region 
NOAA 0960 



the loop top turns out to be ~ 1.5 times above that 
of the corresponding field line. This field line (the 
first entry in table[2]) again shows a relatively small 
value \a\. Since this a value gave the best fit of 
linear force-free field lines to the loop projection in 
the images, we conclude that the linear force-free 



assumption is often not adequate (cf. Wiegelmann 
etaTpOOll ). 



Discussion and outlook 

We demonstrated that EUV data from the new 
STEREO spacecraft allows for the first time to 
make a reliable stereoscopic reconstruction of the 
spatial distribution of hot, magnetically confined 
coronal plasma and, by inference, provide a full 
three dimensional view of the arrangement of coro- 
nal field lines. We found that linear force-free field 
models are helpful to establish correspondences 
between the loops observed in the STEREO im- 
age pairs. The field lines from these linear force- 
free models need not be physical but only serve 
as a first order approximation to the final loops. 
Realistic magnetic field models of the corona will 
have to be judged by their capability to yield field 
lines in agreement with the stereoscopically recon- 
structed loops. Our scheme to determine corre- 
spondences will become even more valuable when 
the stereo base angle grows and loop structures be- 
come more difficult to be identified in the image 
pairs. 

The reconstructions will also allow more precise 
analyses of emissions from loops. The observed 
brightness of EUV loops is, e.g., strongly modi- 
fied by the inverse cosine of the angle between the 
line of sight and the loop's local tangent. This 
may, besides other effects, contribute to the en- 
hanced EUV brightness of the lower loop segments 
commonly observed on the solar disk: these loop 
segments close to the loop's foot points are more 
aligned with the radial direction and they make 
a small angle with the view direction. This may 
cause them to appear brighter than the loop top 
which is viewed at more or less right angles. 



Other applications have been proposed (As 



chwanden 2005 Aschwanden et al. 20061. E.g., 



the amount of twist of a reconstructed loop in- 
dicates how close the flux tube is to a kink in- 
stability. Torok et al. (20041 found a threshold 



of about 3.57T in numerical simulations for the 
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twist <f> = LB^/rBu. Here L is the length of the 
flux tube, -By the toroidal field along its axis and 
B^ ~ aB\\r/2 the poloidal field at a radius r from 
the flux tube centre. In some cases it may be pos- 
sible to resolve the number of turns n which a field 
line makes about the flux tube centre from stereo- 
scopic reconstruction and thus to determine the 
twist from $ = 2ixn. Likewise, the twist is also 
related to a and L by $ = a.L/2. For the ac- 
tive region observed here, table [2] gives values of 
<f> < 0.5 well below the kink instability threshold. 

Another perspective for stereoscopic loop re- 
construction is the analysis of loop oscillations 
from a series of image pairs. The reconstructed 
loops will allow us to determine the transverse po- 



larisation of these oscillations (Aschwandcn ct al 



2002 Wang fc Solanki 20041. Since the coronal 



magnetic field has a complicated geometry with- 
out symmetries, the frequency of these oscillations 
will significantly depend on this polarisation. Note 
that these phenomena are invisible in the mag- 
netic surface data and therefore cannot be re- 
trieved from field extrapolations, which in addi- 
tion require a stationary magnetic field. 
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